MATH 18.152 COURSE NOTES - CLASS MEETING # 10 



18.152 Introduction to PDEs, Fall 2011 Professor: Jared Speck 

Class Meeting # 10: Introduction to the Wave Equation 



I. What is the wave equation? 
The standard wave equation for a function u(t,x) (where t e M., x e M. n ) is 



(1.0.1) -\d?u + Au 



(1.0.1) is second order and linear. The constant c > is called the speed (this terminology will be 
justified as our course progresses), and it has dimensions of ~fi^r • Note that heuristically speaking, 



if we let c -»• oo, then (1.0.1 ) becomes Laplace's equation. However, as we will see, in order to have 



a well-posed problem for (1.0.1), we will need to specify Cauchy (i.e. initial data) for u and also 
d t u. The fact that we need to specify Cauchy data is in stark contrast to Laplace's equation, but 
is analogous to the heat equation. The fact that we need to specify two pieces of Cauchy data is 
connected to the fact that the wave equation is second order in time. 

2. Where does it come from? 



Equation (1.0.1) arises in an incredible variety of physical contexts, especially those involving 
disturbances that propagate at a finite speed. Let's discuss how the wave equation arises as an 
approximation to the equations of fluid mechanics. For simplicity, let's only discuss the case of 1 
spatial dimension. The equations of fluid mechanics, which are known as the Euler equations, take 
the following form in 1 + 1 dimensions: 



(2.0.2a) 
(2.0.2b) 



d t p + d x (pv) =0, 
d t (pv) + d x (pv 2 ) = -d x p, 



where p{t,x) is the fluid mass density, v(t,x) is the fluid velocity, and p(t,x) is the pressure. 
Equation (2.0.2a) implies the conservation of mass, and equation (2.0.2b) is Newton's second law: 



the rate of change of fluid momentum is equal to the force, which is created by the pressure gradient 
(i.e., the -d x p term). The Euler equations are highly nonlinear, and we are very far from obtaining 
a full understanding of how their solutions behave in general. 

A fundamental aspect of fluid mechanics is that the system is not closed because there are not 
enough equations. A common method of achieving closure is by choosing an equation of state, which 
is a relationship between the fluid variables. This relationship is often empirically determined. A 
commonly studied equation of state is 



(2.0.3a) 



p = Kpt 



where 7 > 1 and K > are constants. For future use, we note that under (2.0.3a), we have 
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(2.0.3b) 
(2.0.3c) 



d x p = K-ip 1 1 d x p, 

d 2 xP = K lp i- l dlp + ^ 7 (7 - l)p"- 2 {d x p)\ 



Also for future use, we differentiate (2.0.2a) with respect to t and (2.0.2b) with respect to x to 
deduce that 



(2.0.4a) 
(2.0.4b) 



d 2 t p + pd t d x v + vd t d x p + d t pd x v + d t vd x p = 0, 
pd t d x v + vd t d x p + d t pd x v + d t vd x p + dip + Apvd^v + 2d x pd x v = -dip. 



The theory of acoustics is based on linearizing (i.e. throwing away the nonlinear terms) the 
equations (2.0.4a) - (2.0.4b) around the static solutions p = p = const > 0, v = 0, p = p = const > 0. 



These static solutions describe a fluid at rest. Let's assume that we make a small perturbation of 
this solution, i.e., that v is small, and that 



(2.0.5) 

where 5(t,x) is a small function. 



p = p + 5, 



Using the expansion (2.0.5), we now throw away (with the help of (2.0.3c)) all of the quadratic and 



higher-order small terms from (2.0.4a) - (2.0.4b) to obtain the following approximating system 



(the quantities that are assumed to be small are v, 5, and all of their partial derivatives): 



(2.0.6a) 
(2.0.6b) 



d 2 t 5 + pd t d x v = 0, 

pd t d x v = -K^p^dlS. 



Comparing (2.0.6a) and (2.0.6b), we see that 5 verifies the following approximating equation 



(2.0.7) 



-d 2 t 5 + Kjp^dld = 0. 



Equation (2.0.7) is a wave equation for the perturbation S(t, x)\ It models the propagation of sound 



waves. This is the linear theory of acoustics! Note that the speed associated to the equation (2.0.7) 
depends on the background density p : 



(2.0.8) c = \lK- i p"<- 1 . 

When 7 > 1, higher background density faster sound speed propagation. 

Remark 2.0.1. For air under "normal" atmospheric conditions, 7 = 1.4 is a pretty good model. 

3. Some Well-Posed Problems 

Recall that well-posed PDEs have three important properties: 

• Given suitable data, a solution exists. 

• The solution is unique. 

• The solution depends continuously on the data. 
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Perhaps the most often studied well-posed problem for the wave equation is the global Cauchy 
problem in 1 + n spacetime dimensions: 



(3.0.9a) 
(3.0.9b) 
(3.0.9c) 



-dfu(t,x) + A x u(t,x) = 0, (t,x)e 
u(0, x) = f(x), x e ] 
d t u(0, x) - g(x), x € I 



We now mention some additional well-posed problems in the case of 1 + 1 dimensions. We assume 
that u verifies the wave equation for (t,x) £ (-00,00) x [0,L] and that Cauchy data is given: 



(3.0.10a) 
(3.0.10b) 
(3.0.10c) 



-dfu(t,x) + dlu(t,x) =0, (t,x)eRx [0,L], 
u(Q,x)=f(x), xe[0,L], 
d t u(0, x) = g(x), xe[0,L]. 



Unlike in the case of (3.0.9a) - (3.0.9c), because of the finiteness of the interval [0, L], we need 
to supplement (3.0.10a) - (3.0.10c) with additional conditions in order to generated a well-posed 



problem. Here are some well-known ways of generating a well-posed problem; they are essentially 
the same as in the case of the heat equation. 

(1) Dirichlet data: also specifying u(t,0) - a(t), u(t,L) = b(t) for t > 

(2) Neumann data: also specifying d x u(t, 0) = a(t), d x u(t,L) = b(t) for t > 

(3) Robin data: also specifying d x u(t, 0) = ku(t,0), d x u(t,L) = -ku(t,L) for t > 0, where k > 
is a constant 

(4) Mixed data: e.g. one kind of data at x = 0, and a different kind at x = L 

4. 1 + 1 SPACETIME DIMENSIONS 

Let's consider the wave equation with speed c in 1 + 1 dimensions: 



(4.0.11) 



-c d*u(r, x) + d x u(r, x) = 0. 



Let's first note the following fact: if /, g are any differentiable functions, then u(x,t) d = f(x-cr) 



and u(x,t) d = g(x + cr) solve (4.0.11). The first is called a right-traveling wave, and the second is 
called a left-traveling wave. To visualized wave propagation in 1 + 1 dimensions, you can imagine 
that the graph of /(■) and g(-) are translated to the right/left at a speed c. This gives a good idea 
of what wave motion looks like in 1 + 1 dimensions. In particular, the amplitudes of the traveling 
wave solutions are preserved in time. As we will see, wave propagation in higher dimensions is quite 
different. In higher dimensions, the amplitudes decay in time due to the spreading out of the waves. 
You will study the case of 1 + 3 spatial dimension in one of your homework exercises; you will show 
that in this case, the amplitudes decay at a rate of order t^ 1 as t -*■ 00. 



Remark 4.0.2. Not all wave solutions in 1 + 1 dimensions are traveling waves; see Theorem 4.1 



By making the change of variables t A = ct, we can transform equation (4.0.11 ) into a wave equation 
with speed equal to 1 : 



(4.0.12) 



-dtu(t, x) + dlu{t, x) = 0. 
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This makes our life a bit easier. Let's now consider the global Cauchy problem by supplementing 



(4.0.12) with the initial data 



(4.0.13) 



u(0,x) = f(x), 



d t u(0,x) = g{x). 



As we will see, (4.0.12) + (4.0.13) has a unique solution that has a nice representation. 



Theorem 4.1 (d'Alembert's formula). Assume that f £ C 2 (1R) and g e C^IR). Then the 
unique solution u{t,x) to (4.0.12) + ( |4.0.13[ ) satisfies u e C 2 ([0, oo) x R) and can be represented by 
d'Alembert's formula: 



(4.0.14) 



u(t,x) 



g(z)dz. 



Remark 4.0.3. For the wave equation -c 2 dfu + d 2 u = formula (4.0.14) is replaced with 
(4.0.15) u(t,x) ' ' 



f(x + ct) + f(x-ct)) + — / g(z)dz. 

' IC Jz=x-ct 



Remark 4.0.4. Equation (4.0.14) illustrates the finite speed of propagation property associated 



to the wave equation. More precisely, the value of the solution at (t,x) is only influenced by the 



"initial data interval" {(0,y) | x - t < y < x + 1}; changes to the initial data (4.0.13) outside of this 



interval have no effect on the solution at (t,x). We will reexamine this property later in the course 
with the help of energy methods. 



Proof. To derive (4.0.14), it is convenient to introduce a change of variables called null coordinates: 



def , 
q = t - x, 



(4.0.16) 
(4.0.17) 

The chain rule implies the following relationships between partial derivatives: 



def , 

s = t + x. 



(4.0.18) d g = -(d t -d x ), d s = -(d t + d x ), 

(4.0.19) d t = d q + d s , d x = d s -d q . 

The operators d q and d s can be viewed as directional derivatives in the (t,x) Cartesian spacetime 
direction .5(1,-1) and .5(1,1) respectively. These null directions, which are sometimes called 
characteristic directions, are extremely important. In the future, we will discuss the notion of a 
characteristic direction in a general setting. 



It is now easy to see that (4.0.12) takes the following form in null coordinates: 



(4.0.20) 



d s d q u = 0. 



Integrating (4.0.20) with respect to s, we have that 



(4.0.21) 

where if is a function of q. 



d q u = H(q), 
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Note that the value of q is the same for the pair of Cartesian spacetime points (r, y) and (0, y-r). 
Thus, using the initial conditions (4.0.13), we have that 



(4.0.22) 



d q u(r,y) = d q u(0,y-r) = - d x )u)(0,y - r) = -(g(y -r) - f'(y -r)). 



Similarly, interchanging the partial derivatives in ( 4.0.20[ ) to deduce d s d q u = 0, we conclude that 



d s u(T,y) = -{g{y + r) + f'iy + r)). 



(4.0.23) 

Adding ( |4.0.22| ) and ( |4.0.23[ ), and using ( |4.0.18| ), we have that 
(4.0.24) d t u(t,x) = Uf{x + 1) - f{x -t) + g{x + t) + g(x - *))• 



Integrating (4.0.24) in time with respect to t from to t, and again using the initial conditions 



(4.0.13), we have that 



/(*) 



(4.0.25) u(t, x) = u(0, x) +^(/ + t) - f{x) + f{x - 1) - /(a;)) + ^ J gix + r) + gix - r) dr 

1 . \ 1 rz=x+t 

= i;[f(x + t) + f(x-t)) + - giz)dz, 

2 v / 2 Jz=x-t 

where to derive the last equality, we made the integration change of variables z = x + r for the 
g(x + t) term, and the change of variables z = x - t for the gix - r) term. We have thus derived 



(4.0.14). 



□ 



Without a lot of additional effort, we can extend Theorem AA_ to apply to the following initial 
+ boundary value PDE in 1 + 1 dimensions; the result is stated and proved in the next corollary. 
This PDE would arise in the study of e.g. the following idealized problem: a description of the 
propagation of waves on an infinitely long vibrating string with one end fixed. Furthermore, the 



corollary will later play a role in our extension of Theorem 4.1 to the case of 1 + 3 dimensions. 



Corollary 4.0.1. Let f e C 2 ([0,oo)), g e (7 1 ([0» oo ))» and assume that /(0) = g(0) = 0. Then the 
unique solution to the following 1 + 1 dimensional initial + boundary value problem 



(4.0.26a) 
(4.0.26b) 
(4.0.26c) 
(4.0.26(1) 



-d^uit, x) + d^uit, x) = 



(t,x)e[0,oo)x(0,oo), 
w(t,0) = 0, te[0,oo), 
u(0,x) = fix), xe(0, oo), 
d t u(0,x) = g(x), xe(0,oo) 



satisfies u e C 2 ([0, oo) x [0, oo)). Furthermore, it can be represented as 



(4.0.27) 



uit,x) 



' l(f(x + t)+fix-t)) + lf*:* + _ t tl giz)dz, if0<t<x, 
^ i(f(x + t)-f(t-x)) + if^g(z)dz, if0<x<t. 
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Proof. The idea is that if we extend u to be odd in x, then we can reduce the problem to the case 



of Theorem 4.1 Motivated by this, we define 



(4.0.28) 


u(t,x) d = < 


r u(t,x), if t > 0,x > 0, 
^ -u(t,-x), if t > 0,x < 0, 


(4.0.29) 




r f(x), ifx>0, 

[ -fi-x), ifx<0, ' 


(4.0.30) 


~x >. def 

g(x) = < 


gix), if x > 0, 
-gi~x), if x < 0. 



Since x) solves (|4.0.26a|), it follows that u(t, x) is a solution to the wave equation (|4.0.12|) for 
it,x) eR 



(4.0.31) 



with initial data m(0,x) = fix), dtuit,x) = gix). Thus, by (4.0.14), we have that 

uit,x) = -(fix + t) + f(x-t)) + - / g(z)dz. 
2 v / 2 Jz=x-t 



The expression (4.0.27) now easily follows from considering (4.0.31) separately in the spacetime 
regions {(t,x) |0 < t < x} and {it,x) |0 < x < t}, and from the definitions (4.0.28) - (4.0.30); note 
that in the case {(t, x) |0 < t < x}, since ~g is odd, the part of the integral from x - 1 to t -x cancels 
and thus the only net contribution comes from the integration interval [\x -t\,x + t]. 

□ 
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